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Abstract 

We prove a shape theorem for a growing set of simple random 
walks on Z*^, known as frog model. The dynamics of this process is 
described as follows: There are active particles, which perform in- 
dependent discrete time SRWs, and sleeping particles, which do not 
move. When a sleeping particle is hit by an active particle, the for- 
mer becomes active as well. Initially, a random number of particles is 
placed into each site. At time all particles are sleeping, except for 
those placed at the origin. We prove that the set of all sites visited by 
active particles, rescaled by the elapsed time, converges to a compact 
convex set. 
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1 Introduction 



In this note we study a discrete time particle system in if- named frog model. 
In this model there are active particles, which move as independent simple 
random walks (SRWs) on Z'', and sleeping particles, which do not move until 
activated. At time zero there is a random number ri{x) of particles at each 
site X of the lattice, where {r]{x)^x G Z'^} are i.i.d., and all the particles are 
sleeping except for those that might be placed at the origin. Those active 
particles start to perform a discrete time SRW. From then on when an active 
particle jumps on a sleeping particle, the latter wakes up and starts jumping 
independently, also performing a SRW. If the origin is initially occupied then 
the number of active particles grows to infinity as active particles jump on 
sites that have not been visited before, awakening the particles that are 
sitting there. Let us underline that the active particles do not interact with 
each other and there is no "one-particle-per-site" rule. 

The frog model can be viewed as a model for describing information 
spreading. The original idea is that every active particle has some informa- 
tion and it shares that information with a sleeping particle at the time the 
former jumps on the latter. Particles that have the information move freely 
helping in the process of spreading information. The model that we deal 
with in this paper is a discrete-time version of that proposed by R. Durrett 
(1996, private communication), who also suggested the term "frog model". 



The first published result on this model is due to Teles, Wormald pT 
where it was referred to as the "egg model". They proved that, starting 
from the one-particle-per-site initial configuration, the origin will be visited 
infinitely often a.s. Popov proved that the last result holds in dimension 
(i > 3 for the initial configuration with a sleeping particle (or "egg" ) at each 
a; 7^ with probability Q;/||a;|p, a being a large positive constant. In Alves 
et al. a modification of the present model was studied from the point of 
view of extinction and survival. The difference of the model of from the 
model of this paper is that in the former active particles may disappear on 
each step. Recently A. Ramirez and V. Sidoravicius communicated to us 
that they are working on a continuous-time analog of this model, and that 
they have proved some results such as shape theorem and convergence to the 
product of Poissons. 

In Alves et al. it was proved that, starting from the one-particle- 
per-site initial configuration, the set of the original positions of all active 
particles, rescaled by the elapsed time, converges to a nonempty compact 
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convex set. In the present paper we generalize the main result of |T| to the 
case of random initial configuration. It turns out that for the case when the 
initial configuration contains empty sites, this generalization is nontrivial. 

Now we define the model in a formal way. Let (N^ , i/) be a probability 
space, where Bi is the product sigma algebra and u is the translation invariant 
product measure determined by the distribution of {i]{x) : x G Z'^}. 

For each tu e N^" let {{Sl,^{u)) x e Z"', 1 < < u{x)} be the inde- 
pendent simple random walks which are executed by the particles in u when 
they are activated. We define S'g^(a;) = x, for all x E Z'^ and 1 < A; < a;(a;). 
Denote by Q'^ the path space of the trajectories of the random walks start- 
ing from the initial configuration u and by P^^ the corresponding path space 
measure. Let P be the measure on i7 = HcjeN^^'^l^ ^ ^'^) obtained by taking 
the base measure on N^'' to be the product measure u and the conditional 
measure P[ ■ | u;] = P^^. For each u G N^"*, let 

t{x, z){uj) = min{?T, : 5'^;,(co') = z for some k,l < k < uj{x)} 

and 

m 

T{x,z){u) = inf I (1.1) 

i=l 

where the infimum is taken over all the finite sequences x = Xo,Xi, . . . ,Xm = 
z. Note that t{x, z){u) = T(x, z){lj) = oo when lj{x) = 0. Note also that for 
d > 3, t{x, z){uj) = oo with positive probability even when uj{x) > 1. 

Let us define the set of sites which were visited by active particles up to 
time n, provided that initially the active particles were only in x. Namely, 

Ciuj) = {yeZ'':Tix,y)iuj)<n}. 

We are mostly concerned with C,n '■= asymptotic behaviour. In 

order to analyze that behaviour, define 

C(^) = {y + (-1/2, l/2r : y G C(^)} C 

and ■■= C- 

The main result of this paper is the following 
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Theorem 1.1 For any dimension d > 1 there is a nonempty convex set A = 
A{d, v) C such that for v-almost all initial configurations uo, conditioned 
on {77(0) > 1}, we have for any < e < 1 

(1-£)^C ^ C {l + e)A 

for all n large enough V^-a.s. 

Note that, ahhough Theorem |1.1| establishes the existence of the asymp- 
totic shape A, it is difficult to identify exactly this shape. Of course, A is 
symmetric and ^ C S), where, denoting ||x||i = la;^"^-*! + ■ ■ ■ + 

D = {x G M"' : ||x||i < 1}. 

Also, note that if the initial configuration is augmented (i.e. some new parti- 
cles are added), then the asymptotic shape (when it exists) augments as well. 
In the paper it was shown that if the initial configuration is constructed 
by adding m particles to each site and m is large enough, then the limiting 
shape A contains some pieces of the boundary of 2) (a "flat edge" result). 
Now we show that if the distribution of rj is heavy-tailed enough, then the 
limiting shape A coincides with 2) (a "full diamond" result). 

Theorem 1.2 Suppose that for some positive 6 < d and for all n large 
enough we have 

P[r/(x) >n]> {\ogn)-^. (1.2) 



Then, Theorem 1.1 is verified with A = T). 



2 Proofs 

Proof of Theorem \1 . i| . 

Step 1. First, we state a few preparatory results, which concern mainly the 
tails of random variables T(-, ■). Let us begin by recalling a technical fact 
from [Q. 

Lemma 2.1 Suppose that uj{x) = 1 for all a; G Z*^. For all d > 1 and 

Xo G Z'^ there exist positive finite constants ai = ai{xo,d) and Pi = Pi{d) 
such that 

P[T(0,Xo) > rn] < aiexp{— m^^} 

for all m. 
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Proof. Here we give only the main ideas of the proof, as it was given in full 
detail in (Theorem 3.2). 

1. The case d > 4. Pick n > WxqW^, where || ■ || is the Euclidean norm, 
and fix some e such that < e < 2(d~2) • Define for 1 < i < [d/2\ the sets 
pj^) := {x e Z'^ : < mV2+e|. 

1.1. Consider the trajectory of the initial particle until time n. With over- 
whelming probability it stays in 'D'f'l all that time and awakens at least 
0{n^^^) sleeping particles. To see why the last claim is true, divide the 
time interval [0,n] into disjoint subintervals of size n^~^. During a fixed 
subinterval, the expected size of the corresponding subrange is of order n^~^ 
(this can be seen by considering all the sites y which are at distance at most 
yri}^ from the position of the initial particle at the beginning of the time 
subinterval, estimating the probability that y is hit and summing over y; 
or just by applying directly the known results about the size of the range 
of SRW, see e.g. Hughes 0, pp. 333, 338). As the size of the subrange is 
certainly not greater than n^"*^, with probability bounded away from it will 
be 0{n^~^) (this follows from the following fact: For any random variable X 
with < X < a a.s. and EX > 6 it is true that P[X > h/2\ > h/{2a)). So, 
with overwhelming probability the size of at least one of the subranges will 
be of order n^~^ . 

1.2. Pick n^~^ particles in which were awakened by the initial particle 
by time n. Divide them into n'^ disjoint groups of size Observe the 
particles of a fixed group during ra^^^'^ time units after their activation. Let C, 
be the number of sites in V^^l \ ^["i visited by the particles of that group 
along the time interval mentioned above. Clearly, C, < n^~^^ x ra^"*"^^ = 
and the direct computation (for each y G \ V^^, using that the 
particles from the group are independent, we compute a lower bound on 
the probability that at least one particle of the group hits y, and then sum 
over y) shows that E( = Oin^), so with probability bounded away from 0, C, 
is of order n^. Considering now all the yf groups and using the independence 
again, we obtain that with overwhelming probability there will be O(n^) 
particles in T>^2e \ '^i"e t)y time n + n^^'^'^ . 

1.3. Now, proceeding in the same spirit, we use those O(n^) particles from 
\ to awaken O(n^) particles in V^^^ \ V^^, and so on, to get finally 

©(nL'^/^J) g^p^jyg particles in I^[rf/2j,£ \ ^^L^/aj-i.e ^ + 0{n^+^^) with 

overwhelming probability. 

1.4. Considering those ©(nL'^/^-l) particles in 2^1^/21 ^ \^|d/2|-i£ waiting 
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^i+2e yj^ji^g Qf iiYae more, one gets that with overwhelming probabihty at 
least one of those particles will hit Xq by time m^^ := n + 0{n^'^'^'^). 

2. The case d = 3. Again, let n > ||a;o||^. Take e < 1/4 and consider the 
n^^'^ particles in ©["J awakened by the initial particle until time n. Now, as 
in the last step of the argument for d > 4, until time := n + 0{n^~^'^^) 
with overwhelming probability at least one of those particles will hit Xq. 

3. The case d = 2. This case is treated analogously to the case d = 3. 
That is, first, dividing the time interval [0,n] into n^/^ subintervals of size 
n^/^ (i.e., taking e = 1/2) and using the fact that the expected size of the 
range of SRW grows like 0{n/ log n), we prove that with large probability the 
original particle will awaken 0(n^/^/logn) sleeping particles in the ball of 
radius n until the moment n, where n > ||xo|p. Considering the independent 
random walks of those particles until time m2 := n + 0(n^), we get the result. 

4. The case d = 1. This was not treated in Theorem 3.2 of |l|], but anyway 
it is quite analogous to the cases d = 2,3. First, by time n we will have n^^^ 
active particles situated not farther than n from the origin. Then, waiting 
until mi := n + n"^ one gets the result. □ 

For d > 3 denote Ed = {Q{d — 2))~^; clearly, when doing the proof of 



Lemma |2.1| in dimension d > 3, one may fix e = Ed- From the above proof 
one can deduce that there exist deterministic constants hd (which depend 
only on dimension) such that rrid < hd'n^^'^^'^, d > 3, rrid < hdU^, d = 1,2. 
This means that, in order to obtain an upper bound on P[T(0, sq) > m], one 
should follow the steps of the proof of Lemma |2.1| with n = (m/hdY^^^'^'^^^ 



d > 3, or n = {m/hdY^'^, d = 1, 2. Keeping this in mind and denoting 

where u is the initial configuration and i? is a finite subset of Z"^, consider 
the following 

Definition 2.1 Let pi = P[?7(0) > 1], andm he any positive integer. A fixed 
initial configuration uj is called m-good if 



d> A and 



for any ball B of radius nd{mY^ ^'^^/'^ which is fully inside X)^"^^"^^^ 
we have 71^{B) > pi/2; 
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- we have 7^.(p£f'")) \ P^S^) > Pi/2 for alH = 2, . . . , [d/2\; 

• d = 3, and for any ball B of radius n^{m)^^~'^^^/'-^ which is fully inside 
Pg'")) we haven^{B)>pj2; 

• d = 1,2, and for any ball B of radius nci{mY^^ situated not farther 
than na{m) from the origin, we have TZ^{B) > pi/2, 

where nd{m) = {m/hciY^^^~^'^^'^\ d > 3, or nd{m) = [m/hdY^'^ , d = 1,2, and 
the notation T>^^^^ is from the proof of Lemma \2. 4 



Lemma 2.2 For all d > 1 and xq G Z'^ there exist positive finite constants 
«2 = (y2{d,pi) and P2 = P2{d,pi) such that ifu is m-good and ||xoP < ndijn) 
(the notation nd{m) is from Definition \2. then 

P^[r(0,Xo) > m] < a2exp{-m^2|^ 
Proof. From Definition ^.1\ it follows that, by following the steps of the proof 



of Lemma |2]1| one can prove Lemma p.2| . For example, in the part 1.1 the 
expected amount of sleeping particles activated during a fixed time subinter- 
val differs from the expected size of the subrange only by a constant factor 
(depending only on pi and d) when uj is m-good. In the part 1.2, will be of 
the same order for m-good uj and for one-particle-per-site uj. In general, each 
time that one is computing the expected number of newly awakened particles 



in the proof of Lemma p.l| , considering m-good initial configuration instead 
of one-particle-per-site initial configuration costs only a constant factor. An- 
other observation, which is crucial for the subsequent discussion, comes after 
examining the proof of Lemma p.l| : As long as ||xo||^ < n, the estimates that 
one gets on P[T(0,Xo) > mj are uniform in Xq. □ 

Define hd{m) to be equal to nd{mY^~'^'^'^^'>f'^ , d > 3, and to nd{my^^, 
d=l,2. 

Lemma 2.3 For all d > 1 there exist positive finite constants = as{d,pi) 
and /?3 = (3^{d,pi) such that 

Y*[r] is m-good \ rj(xi) = ■ ■ ■ = ri(xko) = 0] > 1 — «3 exp{— m^^} 

for any fixed collection of sites Xi, . . . , Xkg G Z'', < hd{mY- 
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Proof. Note that 'ha{m)'^ is small relative to the sizes of sets which were 
considered in Definition |2.1| . As r){x) are i.i.d., it is straightforward to prove 
this fact by using the large deviations technique. □ 

In the sequel we will need the following basic fact which is stated without 
proof here. 

Lemma 2.4 Let {Xi,i > 1) be a sequence of nonnegative random variables, 
and there exist (pn, n = 1,2,..., such that supj P[Xj > n] < and J2n Vn < 
oo. Then X^jn — > a.s. 

Step 2. Now the goal is to verify the conditions of Kingman-Liggett subad- 
ditive ergodic theorem 0, § for the sequence of random variables y (m, n) 
defined below. 

For a fixed x G Z'^ and an e satisfying the condition co'(O) > 1 
define a sequence of positive integers {f^'j^Q as follows: 

vl = 0, 
vl_^i = min{?T, > : u{nx) > 1}. 

In words, v^x, k = 0,1,2,..., are the sites on the ray {nx,n > 0} 
which are occupied in the initial configuration. Clearly, for any x G Z*^ 
and z = 1, 2, . . . it holds that P[ff — = fc] = pi(l — pi)^~^ (pi is from 



Definition 2.1). Now, for m,n > consider the collection of random variables 



Y{m,n) =T{v^x,v^x). 

It is important to observe that the random variables {T{x,y) : x,y E Z'^} 
are subadditive in the sense that for all u, 

T{x, z) < T{x, y) + T{y, z) for all x,y,ze . (2.1) 

Here is the explanation. Fix the initial configuration uj. If the site y was 
empty in the initial configuration, then T(|/, 2;) = oo and ( |2.1| ) follows. Now, 
suppose that ojiy) > 1. If site z is reached before site y, then (pTT| ) is evident. 
If that does not happen, recall that the random variables T{y, z),y, z E 'Z^'^ are 
constructed using the same collection of the random variables S^j^, i.e., each 
particle follows the same trajectory as soon as it wakes up. So the process 
departuring from only site y awakened (the one which gives the passage time 
T{y,z)) is coupled with the original process (i.e., that started from x), and 
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for the latter one may have other particles awakened at time T(x, y) besides 
that from y. Consequently, from (|1 . 1|) it follows that T(x, z) — T{x, y), which 
is the remaining time to reach site z for the original process, is less than or 
equal to T{y,z), thus proving (pTTp. The equation ( |2.1| ) shows that for all 
initial configurations u, 

Y{m, k) < Y{m, n) + F(n, k) for all m,n,k e N. (2.2) 

Now, let us verify the conditions of Kingman-Liggett subadditive ergodic 
theorem for the random variables Y{m,n). 

Since i/ is a product measure we can easily establish the following result: 
Given nonnegative integer numbers mi, m2, . . . , m„, p, the joint distributions 
of random variables {¥{7711,7712), .. . ,y(m„_i,m„)} and {y(mi +p, m2 + 
p), . ■ ■ , y(m„_i +p, m„ +p)} are equal. From this we obtain the stationarity 
of the sequence {¥{{71 — l)k,7ik)}^^^ for each k E N. Ergodicity of the 
sequence follows from the observation that the events {T(f ^^^.x, v^^^_^_^'^i^x) < 
a} and {T{v^^i^x,v^^^_^_^.^f^x) < b} are independent provided a + b < \\{ni — 
n2)kx\\i. 



Next we prove that ET(0,ffx) < 00. Using Lemmas |2.2| and |2.3| , for y 
such that \\y\\ < hd{7nY^'^ (note that ^^(m) < 7id{7n)) and for any collection 
of sites -21,... ,ZkQ G Z'^, ko < fidimY, we have, denoting P*[-] = P[ ■ | 
r^{zi) = --- = 7]{zko) = 0], 

P*[T{0,y) > 7n] < P*[T{0,y) > m | 77 is m-good] + P*[7] is not m-good] 
< Q;2exp{— + asexpj— m'^^}. (2.3) 



Using ( |2.3| ), we get 



P[T(0, <a;) >m] = ^ P[< = k]P[T{0, kx) > m \ = k] 

k=l 

< P[< = fc]P[T(0, fcx) > m K = k] 

k<\\x\\~^hj(m)^/'^ 

+ E pi(i-pi)'"' 

fc>||x||~-'-n£j(m)-^/^ 

< a2 exp{— m^^} + 0:3 exp{— m^^} 

^ (^i _ p^yM-'Mmy/\ (2.4) 

so Er(0, 1) = J2m>i P[^(0, <a;) > m] < 00. 
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Thus, we have verified the conditions of Kingman-Liggett subadditive 
ergodic theorem for the sequence {Y{m,n) = T{v^x,v^x)}. Therefore, one 
gets that there exists fj,'{x) such that 

hm ^ ' ^ = fi'{x) < oo P-a.s. (2.5) 

n—*oo fl 

Step 3. Now the goal is to pass from the random sequence {vfx, > 0) to 
the whole ray {kx, k > 0). 

Denote := Let n G N satisfy < n < v^f^^y Then, by 

the subadditivity, 

T(0,nx) ^ T(0,t;g(„)X)A;(n) , T(t;g(„)X, nx) 
n k[n) n n 

Note that, since is a product measure, f^^-^-j — n has the same distribution 
as vf. This means that the upper bound on P[T(t;^(^^^a:, nx) > m] is also 
given by ( p^.4| ), so by Lemma 2]^ we have that n~^T{v^^^-^x,nx) 0, P-a.s. 



Together with (|2.5| ) and the fact that k{n)/n — > pi P-a.s., this shows that 
one gets from ( ^.61) that 

limsup^^^^^^ <^(x) P-a.s. (2.7) 

n— >oo 

Now, let us prove that 

liminf ^^^'""^^ > /i(x) P-a.s. (2.8) 

For a fixed site y G Z"^, y ^ 0, let [/^ be a random variable defined in the 
following way. At the moment T(0, y) consider the active particle situated 
in y (if there are several such particles, choose one of them by randomization). 
Let Zl, k > T{0,y) be the subsequent walk of this particle. We define r = 
min{fc : //(Z^) > 1} and so Uy := t — T{0,y) is the time that the particle 
needs to travel from y to some site that initially contained sleeping particles. 
Clearly, Uy = corresponds to the case ri{y) > 1. Now, let (n = ^t{o nx)+u„^'^ 
note that i]{(n) > 1- By subadditivity one can write 

T(0, <(„)X) < T(0, nx) + Un. + T(C„, <(„)x). (2.9) 
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To proceed, we need to get some good bounds on the tails of Unx and 
T{(n,vl(^n)^)- Note that there exists a positive constant P which is not de- 
pendent on the dimension, such that P[|{Zj'^',0 < j < k^^"^}] > k^^^] > /3, 
for all > 1. This shows that 



P[|{Z"^ < J < A;}| > A;^/^] > 1 - (1 - /5) 



fcl/2 



So, as any site in {^""^,0 < j < k} initially contained a sleeping particle 
with probability pi, it holds that 

P[Un. >k]<{l- (3f + (1 - p,f'\ (2.10) 
Now, Cn|| < Un, (12.101) implies that 

P[\Kn)^-L\\>k] < P[||<(„)X-nx|| > A;/2]+P[||nx-Cn|| > A:/2] 

< (1 _ + (1 _ ^)(fc/2)V^ ^ (1 _ p^)(fc/2)V* 

=: V^(A;), 

and, proceeding similarly to ( |2.4| ), we get that 



P[T(4,^^fc(„)a;) > m] < a2exp{-m^'} + a3exp{-m^=^} + 7/'(||a;|| ^rid(m)^/^). 

(2.11) 

Thus, by ( CT) and (|2TT|) , dividing (|2;|) by n and applying (|2j) and 
Lemma p.4| , we get ( p.8|) . Now, from (|2.7| ) and ( p.8|) we finally obtain that, 
for all X G Z*^, there exists fi{x) > such that 

lim ^l^^i!^ ^ ^(a;) p.a.s. (2.12) 

Notice that the equation ( p. 121) is already enough to get the proof of 
Theorem |TTT] in dimension 1 with A = [— (/i(l))~^, (yu(l))^^]. So, from now 
on we shall concentrate on the case d > 2. 



Step 4. To proceed with the proof of Theorem |1 . 1| , one has to assure that 
grows at least linearly. For y & M.'^ and a > denote D{y,a) = {x G M*^ : 
\\x - y\\i < a}. 

Lemma 2.5 For all d > 2 there exist constants < 5 < 1, 0:4 > 0, > 0, 
which depend only on the dimension, such that 

'P[D{x,n6) C f„+T(o,x)] > 1 - a4exp{-?2'^''} 

for all n and x, conditioned on {?7(0) > 1}- 
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Proof. The proof of this fact follows the spirit of the proof of Lemmas 4.2 
and 4.3 of The difficulty that arises here is that, for u such that uj{x) = 0, 
T(x,y){uj) = oo for any y, so a direct application of the method of [Q] does 
not work. To get around that difficulty, we use the construction similar to 
that used in Step 3. Recall the definition of random variable Uy from Step 3, 
and define for y ^ 0, T*{y,z) = Uy + T{(y,z), where, similarly to Step 3, 
is the site with ^^(C^) > 1 at which the active particle which was in y at 
time T{0,y) arrived after Uy steps. Now, let y be such that \\x — y\\i = n, 
and let x = yo,yi, . . . ,yn = y he a path connecting x to y such that for 
all i, \\yi — yi-i\\i = 1; note that \\x — y^Wi = k, k = 0, ... ,n. Denote 
Yi = T*{yi_i,yi). From the above construction it follows that 

r,,l/2 



where 



TiO,y)-TiO,x)<J2Y^ = J2' 



j:i+jn^^'^<n 



1, . . . ,?T,^/^. Analogously to ( |2.10|) and (|2.11|) , it is not difficult to get 



that for some constants a^, /?5 > and for all m 

P[Yi >m]< a5exp(-m'^5), (2.13) 

uniformly in x, y and paths connecting them. Consider the event B = {Yi < 
n^/^/2, i = 1, . . . , n}; clearly, from ( p.l3|) it follows that for some a^, (3q > 0, 



P[B] > 1 — ae exp(— n^"^). Note that, if \i — j\ > n^^^, then the random 
variables Y^Ib and YjIb are independent, since if the event B occurs, then 
the random variables Y^ and Yj depend on disjoint sets of random walks. 
So, when B occurs, each cij is a sum of independent random variables. Al- 
though Yi are not identically distributed and we cannot guarantee the exis- 
tence of moment generating function for Yi, the condition ( |2.13| ) allows us. 



by applying Theorem 1.1 of Nagaev 0, to obtain that there exists 5o > 
such that cTj < n^/^/^o with probability at least sub exponentially high. This 
shows that, with overwhelming probability, Yl^=i^i — ^/^o- Thus, if y is 
at distance n from x, with overwhelming probability by time T(0,x) +n/6o 
it will be visited. Now, ii d > 2 and < — < n, then there exists 
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z E such that \\z — a;||i = n and H^; — y\\i equals n or n + 1. Since 
r*(x, y) < T*{x, z) + T*{z, y), the resuh follows with 6 = 5q/2. □ 

Step 5. The next step is to prove that /i(-) can be extended to in such a 
way that ;U is a norm in W^. Let us extend the definition of T{x,y) to the 
whole W'- X M'^ by defining 

T{x,y) = mm{n : ?/ G 

where xq G Z'^ is such that x G (—1/2, 1/2^ + xq. From the fact T(0, nx) > 
n||x||i it follows that fi{x) > \\x\\i for all x G Z*^. 

Lemma 2.6 For any a G M+, x, ?/ G Z*^ we /iatJe 

fi{ax) := lim — ^'^ = a/i(x) P-a.s., (2-14) 

and 

li{x + y) < iJ,{x) + n{y). (2.15) 

Proof. First, the proof of ( |2.14| ) basically repeats what was done on Step 3, 
so we omit it. Let us turn to the proof of (|2.15| ). Clearly, instead of ( p.l5| ) it 
is sufficient to prove that for any x,y G Z*^ 

^^'{x + y)<^^'{x)+^'{y) (2.16) 



with /i'(-) defined by (|2.5| ). 

For fixed n > 1 define a random sequence Sk^ni A; = 0, 1, 2, . . . , in the 
following way: 

SO,n = 0, 

Sk+i,n = min{m > Sk,n ■ vi'^V + K^) > 1}^ 
and let z^^^i = Sk^nV + v^x. By the subadditivity, 

T(0, <+^(a; + y))< T(0, <a;) + T«a;, Zn,n) + T{zn,n, vl+y{x + y)). (2.17) 
We have 

n 

<+n^+y)-z„,„ = xY,{{vr' - vt^) - {v^, - vu)) 

i=l 

n 
i=l 
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Fix arbitrary e > 0. Now, (ff+^ - v^^^), {v^ - (si,„ - Si_i^„), i = 

1, . . . ,n, are i.i.d. random variables, geometrically distributed with param- 
eter (1 — pi). So, in the right-hand side of the above display we have two 
sums of i.i.d. random variables satisfying the Cramer condition and with 
zero mean. So, by the Large Deviation Principle, with exponentially high 
probability + — < £:(||a;|| + ||?/||)ra. By Lemma with subex- 

ponentially high probability, T(z„^„, + y)) < e6~^{\\x\\ + \\y\\)n (the 

site Zn,n is random, so, to see that Lemma is applicable here, note that 
T{zn,n, Vn'^^{x + y)) has the same distribution as T(0, u), where m is a random 
site constructed as follows: First, let u' be the n-th occupied site on the ray 
along {—y), then, u" is the n-th occupied site on the ray along (— x) beginning 
in u', and finally, u is the n-th occupied site on the ray along x + y begin- 
ning in u"). Note that T{v^x, Zn,n) has the same distribution as T{0,v^y), 
so n~^T{v^x, Zn,n) — > /^'(z/) at least in probability. Dividing ( ^.17| ) by n and 
taking the limit in probability, we get 

fi\x + y)< + fi'iy) + e6'W\x\\ + 

when e 0, we obtain ( p. 161 ). □ 

Now, let us show that fi{x) can be extended to the whole M'' in such a way 
that ( [^.12D still holds. This is done in a standard way: First, if X eQ'^ (i.e.. 



all the coordinates of x are rational), we define /i(x) := jj,{mx)/m, where m is 
the smallest positive integer such that mx e Z'^ (the fact that (P?!^ ) still holds 
for X G Q"' follows from (|2.14|) ). Then, by using the fact that is a norm in ''^'^ 



(this follows from Lemma |2.6| ), it is extended to M.'^. To prove that (|2.12|) 
holds for any x G M^, one can proceed as follows: By the subadditivity and 



Lemma if y G Q*^, then lim sup^^^ n ^T{0,nx) < fi{y) + 6 ^||a; — ?/||, and 



lim inf n^oo ''^ ^T{0,nx) > ^{y) — 6 -^{{x — y\\. Approximating x by vectors 
with rational coordinates, one gets the result. 

Step 6. Now everything is ready to finish the proof of Theorem |1.1| . 

Let A := {x E M.'^ : fi{x) < !}■ The following argument is standard (see 
e.g. [^, |, 0), we keep it to preserve the self-containedness of the paper. 

Denote e' = {1 — e)^^ — 1, and e" = 1 — {1 + s)^^. To prove Theorem |1.1| , 
it is enough to prove that nA C ^(^i^s')n and C nA for all n large 

enough, P-a.s. (and so P^^-a.s. for //-almost all uj). 

Choose a finite set F := {xi, . . . ,Xk} C A such that fi{xi) < 1 for 
i = 1,... ,k, and (with 6 from Lemma p.5|) A C U^^^D{xi, e'6). Notice 
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that ( |2.12| ) implies that nF C for all n large enough, P-a.s. Now, 
Lemma |2.5| and Borel-Cantelli imply that P-a.s. for all n large enough we 
have D{nxi, ne'5) C ^^^/, for alH = 1, 2, . . . , fc. So nA C and this 

part of the proof is done. 

Now, let us choose G := {j/i, . . . ,yk} C 2A \ ^ in such a way that 
2A\A C \J^^iD{yi, e"6). Analogously, we get that nGCi^n = for all n large 
enough P-a.s., and that for all n large enough, if fl (2^\^) ^ 0, then 

n nG 7^ 0. This shows that C nA for all n large enough, P-a.s., 



and so concludes the proof of Theorem □ 
Proof of Theorem \1.2^ . Denote 

'Dn = {x eZ'^ : \\x\\i < n}. 

Choose 6 < 1 such that 5 < 6d. Start the process and wait until the mo- 



ment n . By Lemma 2]5| there exist Ci, 70, such that with probability at 
least 1 — expjn'''"} all the frogs which were initially in the ball of radius Gin^ 
centered in will be awake. Clearly, the inequality ( |1.2| ) implies that 

P[r/(a;) < (4d)"] < 1 - {\ogAd)-^n-^ . 

As the number of particles in that ball is of order n^*^, one gets that with 
probability at least 1 — exp{— C2n^'^~''} at time one will have at least one 
activated site x with ri{x) > (Ad)"' and ||x|| < Cin^. Note the following simple 
fact: If X contains at least {2d)"' active particles and — < n, then until 
time n with probability bounded away from at least one of those particles 
will hit y. Using this fact, as x really contains at least 2" groups of (2c?)"' 
particles, we get that with overwhelming probability all the particles in the 
diamond S)n-n» will be awake at time + n, which completes the proof of 



Theorem 1.2. □ 



3 Remarks about continuous time 

A continuous-time version of the frog model can also be considered. Here we 
would like to remark that in the continuous-time context and for the case of 
bounded rj, Theorem |1.1| also holds and its proof can be obtained by following 
the steps of our proof for the discrete case. The difficulty that arises is that, 
for continuous time, it is not evident that (defined by (|2.12|) ) is strictly 
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positive for x 7^ 0, i.e., we must rule out the possibility that the continuous- 
time frog model grows faster than linearly. To overcome that difficulty, note 
the following fact (compare with Lemma 9 on page 16 of Chapter 1 of [^]): 
there exist a positive number (3 such that, being ||a;||i > /3n, P[T(0,a;) < n] 
is exponentially small in n. This fact in turn follows from a domination of 
the frog model by branching random walk. So, we conclude that for a.s. 
bounded t] our method works well in the continuous-time context too. 
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